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ABSTRACT 
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analysis  and  design  of  time  invariant  unity  feedback^  control 


systems  is  presented,  using  cascade  or  feedback  or  both 
types  of  compensation. 

By  using  this  program,  the  user  is  freed  from  the 
tedious,  time  consuming  and  error  prone  method  of  hand 
calculations,  letting  the  computer  handle  these  tasks 
efficiently  and  speedily.  The  user  can  then  concentrate 
fully  on  the  placement  of  poles  and  zeroes  of  the 
compensator (s)  used. 

Design  of  control  systems  by  classical  methods  being 
essentially  a  repetitive,  trial  and  error  procedure,  this 


program  greatly  cuts  down  the  turn  around  time  and  leads  to 
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I. 


INTRODUCTION 


During  the  past  two  decades,  the  scientific  and 
engineering  communities  have  witnessed  an  ever  increasing 
role  of  the  digital  computer  in  the  fields  of  research, 
development  and  analysis  of  systems.  The  computer,  today, 
is  being  used  to  solve  engineering  problems,  whose  solution, 
until  very  recently  required  long  and  tedious  procedures. 
Still,  it  is  probably  true  that  this  machine  has  potential 
not  fully  recognized  yet,  which  is  why  the  attention  of  so 
many  computer  scientists/engineers  is  focused  on  devising 
more  efficient  and  innovative  operating  procedures. 

Control  system  design  is  one  area  where  classical 
theory  has  been  extensively  developed  and  used.  It  is  fair 
to  say  that  even  today,  most  analysis  and  design  problems  of 
linear,  time  invariant  control  systems  can  still  be 
approached  using  the  methods  developed  by  Bode,  Nyquist  and 
others . 

A  totally  new  approach  to  the  design  of  control  systems 
became  available  with  the  development  of  optimal  control 
theory  and  the  state  variable  analysis.  These  methods  have 
been  intensively  developed  in  the  last  10-15  years,  but  now 
their  weaknesses  have  been  exposed  too.  The  "states"  of  the 
plant  may  not  necessarily  represent  physically  measurable 
quantities,  and  consequently  it  may  not  be  possible  to 


implement  the  results  at  all.  Luenberger's  observers, 
designed  to  overcome  this  problem,  can  at  best  provide 
estimates  of  the  state  trajectory.  Furthermore,  the  optimal 
control  approach  to  design  relies  very  heavily  on  mathe¬ 
matical  manipulation,  providing  little  insight  to  the  actual 
working  of  the  plant;  the  only  input  of  the  designer  being 
the  form  of  the  cost  function. 

An  intelligent  use  of  the  speed  and  information 
processing  ability  of  the  digital  computer,  coupled  with  the 
reliable  features  of  classical  theory  appear  to  be  the  best 
solution  to  the  problem  at  hand.  The  classical  approach  to 
design,  being  essentially  a  trial  and  error  method,  if  the 
order  of  the  system  is  fairly  high,  the  number  of  repetitive 
calculations  and  the  time  required  to  perform  these  calcu¬ 
lations  becomes  prohibitively  large,  the  assistance  of  the 
computer  in  such  problems  becomes  indispensable. 

The  work  in  this  thesis  was  to  develop  an  interactive, 
user  oriented  computer  program  that  would  prompt  the  user 
to  input  the  transfer  function  and  cascade/feedback  compen¬ 
sators.  The  program  would  then  display  on  the  IBM  3277  - 
Tektronix  618  dial  screen  terminals  the  Bode  Plot  of  both 
magnitude  and  phase.  The  program  could  be  repeatedly  used, 
with  the  user  having  the  option  to  change/modify  the  compen¬ 
sators,  each  time  viewing  the  effect  of  his  modifications  on 
the  screen  until  he  arrives  at  a  satisfactory  solution. 
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CONTROL  ENGINEERING  ANALYSIS 


II . 

A.  GENERAL 

A  continuous  time  control  system  may  be  represented  in 
one  of  the  following  forms: 

a.  Transfer  functions 

b.  State  equations 

c.  System  block  diagrams  or  signal  flow  graphs 
Algorithms  exist  in  almost  every  undergraduate  control 

engineering  text  to  convert  the  system  representation  from 
one  given  form  to  another.  Gianniotis  (Ref.  1)  describes  a 
simple  method  of  converting  from  transfer  function  to  state 
variable  form  in  Chapter  II. 

The  transfer  function  representation  in  its  most  general 
form  is: 


_  m  a  m- 1  .  Q  m-  2  a  c  c 

m  m- 1  m- 2  o 


B  Sn  +  B  ,Sn  1  +  B  _,Sn  2  +  .  B  SC 

n  n- 1  n-2  o 


Usually  the  mathematical  description  of  the  system  is 
found  in  the  transfer  function  form  in  the  literature. 
Analysis  and  design  of  control  systems  by  classical  methods 
also  requires  the  representation  of  the  system  in  this  form. 


Table  1.  (Contd.) 


0 . 464 1 57£  +  Cl 
0  .ai4  1 74E  +  01 
0  .  :>o9  5  79  E  +  Gi 
0  .  o3u  955  E  +  Ca 
U  .o45  945  £+01 
0 . 774  2t»lE  +  01 
0  .  6 5 7 o 92 u+Ui 
0.950  lioE+Oi 
U.10525QE  +  02 
0  .115  591  E+02 
0.  129  155  E  +  02 
0  .  j.4a Q7..E+02 
0 .158  +o9  E  >02 
0 . 1 75  5o7E  >  02 
0  ..9++o5E  +02 
0.215443E+02 
0  .238658E  +  02 
0 .264  3  75  £  +  02 
0  .2  92  863  E  +  02 
0  .324 +21 £+02 
0  .359  3aQt  +  02 
0 .a98 lOoE  +  02 
0  .■+41  0U4E  +  02 
0 .46a  526  E  +02 
0 .541 167E+U2 
0 .5 9*4o3 E +02 
0 .o64 OoOE  +  02 
0 .7.35  6a9  E  +02 
0  .8 1 4  4 10.E  +  02 
0.902  722E  +  02 
0  .999  996 c+02 


-0.  4oo01aE+01 
-O.o  362a3E  +  01 
-0 . 7  oa  934E  +0  i 
-Q.  94+9oo:  +  0i. 
-0.1 1G539E  +02 
-0.  l27iloc  +  U2 
-0.  1443-i.aE  +  0<i 
-0. 1 o2 1 9oE+02 

-  0.  180814E  +  02 
-0. 20Ui95E+02 
-0.  220a+uE+02 

-  0. 2+123 oE  +  Q2 
-0. 262a5oE+u2 

-  U.  285134E+02 
-0.  d080a3E+02 
-0.3314ooE+02 

-  0. 35543  IE  +  02 
-U.a79o07E  +  02 
-0. 4G4556E+02 
-0.  4^9o2-3E  +  G2 
—0. 45+9o2E+G2 
-0.  48052 9E  +  02 
-0.  50O288E  +  02 
-0. 532206E+02 

-  0.  55a2oUc  +  02 
-0 . 5o4423E+02 
-0.  610677E  +  02 

-  0.  6a7 0Q5E  +  O2 
-0. 6oa395E+02 
-0.  o8983 5E  +  02 
-0.  71O317E  +  02 


—  u.luOoudE  +  Ga 

—  0. i oa+  7+E+  U3 
-0.  io7a7oE+0a 
-0.x714a2E+ua 
-0.17  554ut+  6a 
-0.1  7  9odaE+03 

—  0  .  i.  8  4+  a  3  E  +  0  5 
-0. la9iaaE+03 
-o.  1 9  3^  3o  E+ua 
-0.1 40  7a  2t  +  0a 
-0. 1  0  aol 2  t  +  u3 
-0.^.0  6a  o9E  +  Oa 
-0.2l300aE+0a 
— U.2I 7469E+0a 
-0.22 1735E+  0a 
-0 .22  57  76E  +  06 
-0. 22  95  74E+  0a 
-0.2a3l21E+0a 
-0.23o4iot+03 
— 0.2394o2E+0a 
-0.242^67t+0a 
— 0.2448+0E+03 
-Q.c47145E+Qa 
— 6 .2  4 9a +4 E  + 03 
-0.^5ia02E+0a 
-o.25au  8aE+  0a 
-0.2a47u0t+03 
-0. ^aolodE+Oa 
-0.25  749  6E  +  Oa 
-0.258703E+03 
-0.259793E+O8 
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Tabular  Output  of  Example  4.1 
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— 0 .15 542  7  £+05 

—  0. A59055E  +  03 
-0. 14240*6+05 
-0. A  453  DA  E  +  05 
-0. 1+o567c+05 
-0 . A  5  A4  9  3E  +  O5 
-0.A54429C+O5 

—  0.  157+-44C  +  05 
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The  velocity  constant  is  not  to  be  decreased.  Design  a 
phase  lead  compensator  which  will  provide  M^w  less  than  1.5. 
b.  Solution 

The  first  step  is  to  get  the  transfer  function 
in  the  required  form: 

5.0 

G  (s )  = - (4.2) 

0 . 2  IS 3  +  1.0S2  +  1.0S 

The  Bode  Plot  of  the  uncompensated  system  drawn 
using  the  program  is  shown  on  figure  2  and  the  tabular  output 
on  table  1.  The  uncompensated  system  has  a  phase  margin  of 
-20  degrees.  To  achieve  a  M  o*r  less  than  1.5,  a  phase 
margin  of  44  degrees  or  more  is  required.  Approximately  64 
degrees  of  positive  phase  shift  are  therefore  needed.  Two 
sections  of  lead  filter  are  introduced  as  given  below: 

(S/3.0  +  1)  (S/1.5  +  1) 

G  =  - 

C  2 
(S/10.0  +  1) 

The  compensated  Bode  Plot  is  given  on  figure 
3  showing  a  phase  margin  of  35  degrees. 

It  may  be  pointed  out  that  the  final  values 
for  the  lead  filter  poles  and  zeroes  were  arrived  at  after 
3  iterations  and  this  design  problem  was  solved  in  less  than 
15  minutes. 
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IV.  PROGRAM  PERFORMANCE  INVESTIGATION 


A.  GENERAL 

The  program  was  tested  by  solving  several  linear  control 
system  design  problems.  Very  satisfactory  results  were 
obtained  in  all  cases  with  remarkable  efficiency  and  speed. 
The  only  necessary  condition  is  proper  problem  formulation. 
This  is  true  for  any  interactive  computer  program.  Once  the 
program  is  used  a  couple  of  times,  the  user  gets  the 
necessary  familiarity  with  its  working. 

The  example  problems  presented  below  are  used  to 
demonstrate  the  performance  and  capabilities  of  the  program. 
The  examples  can  also  help  the  user  in  formulating  his  own 
particular  problem.  The  examples  were  selected  from  Thaler's 
'Design  of  Feedback  Systems'  {R?f.  2) 

B.  EXAMPLE  PROBLEMS 

1 .  Example  4.1  :  A  Phase  Lead  Network 

a.  Problem  Statement 

A  positioning  system  is  single  loop  with  unity 
feedback  and  forward  transfer  function 

5.0 

G  (s )  =  - 

s (0 . 7s  +  1)  (0.3s  +  1) 
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TITLES  inputs  the  two  lines  of  text  as  headings  for 


the  Bode  Plot. 

In  addition,  the  main  program  handles  the  tasks  of 
displaying  the  tabular  output  and  the  Bode  Plot  of  the  open 
and  closed  loop  response  of  the  system  as  required. 
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C.  PROGRAM  DESCRIPTION 

The  main  program  is  the  coordination  center  which 
controls  the  calling  of  the  supporting  subroutines,  in  order 
to  input  the  transfer  function,  cascade/feedback  compensators 
and  other  necessary  information.  The  tabular  results  and  the 
Bode  Plot  are  then  displayed  on  the  two  screens  respectively. 

The  main  program  as  well  as  the  accompanying  subroutines 
can  be  found  in  Appendix  A.  They  contain  a  sufficient  number 
of  comment  cards  to  be  self  explanatory. 

The  following  is  a  brief  description  of  the  performance 
and  purpose  of  the  various  subroutines. 

NUMER  inputs  the  numerator  of  the  plant  transfer 
function. 

DENOM  inputs  the  denominator  of  the  plant  transfer 
function. 

CAS CAD  inputs  up  to  6  first  order  lead/lag  filters. 

SECAS  inputs  the  numerator  of  the  second  order  band 
pass/band  stop  filter. 

SECASD  inputs  the  denominator  of  the  second  order 
filter . 

FETCH  determines  the  value  of  the  radial  frequency, 
w,  at  the  origin  of  the  x-axis. 

DECADE  determines  the  number  of  decades  of  frequency 
to  be  spanned. 

FEEDBK  inputs  the  various  parameters  of  feedback 


compensators . 


In  brief,  the  whole  BODPLT  package  works  as  follows: 

A  user,  logged  into  the  VM/CMS  environment  of  the  system  from 
the  dual  screen  terminals,  issues  the  command  DISSPLA  BODPLT. 
The  package  then  assumes  control.  The  program  begins  its 
execution  by  interrogating  the  user  and  calling  the  appro¬ 
priate  subroutine  accordingly.  All  programmer-composed 
subprograms  are  included  in  the  main  program  titled  BODPLT 
FORTRAN . 

The  BODPLT  program  has  the  following  important  features: 
runs  of  the  VM/CMS  time  sharing  system, 
interrogates  the  user  in  entering  all  problem 
specifications  from  the  terminal. 

can  handle  up  to  a  ninth  order  plant  transfer 
function,  six  first  order  and  one  second  order  cascade 
filters. 

prompts  the  user  to  input  the  parameters  of 
velocity  feedback,  acceleration  feedback  or  approximate 
acceleration  feedback  as  required. 

provides  the  solution  in  tabular  form  on  the  IBM 
3277  screen  and  the  Bode  Plot  on  the  TEK  618  terminal. 

can  provide  hard  copy  version  of  the  problem 
specifications  and  tabular  output  by  using  the  RECORD  ON/ 
RECORD  OFF  execs,  and  of  the  BODE  PLOT  on  the  Tektronix 
printer  where  installed. 

allows  problem  specifications  to  be  changed  between 

runs . 


III.  PROGRAMMING  CONSIDERATIONS 


A.  GENERAL 

Any  interactive  software  package,  such  as  the  one 
developed  in  this  thesis  must  provide  a  simple  yet  unambigous 
means  of  data  input  and  output.  The  data  must  be  easy  to 
interpret  and  apply  to  the  problem  at  hand.  Programs 
producing  highly  satisfactory  results  but  requiring  long 
studying  time  and/or  special  programming  skills  are  of 
limited  use  only. 

The  intent  of  this  thesis  is  to  present  such  a  program. 
Computer  Aided  Design  of  Linear  Systems.  Special  care  has 
been  taken  to  develop  the  program  so  that  the  user  has  to 
invest  very  little  time  learning  to  use  it. 

B.  MAIN  FEATURES  OF  THE  PROGRAM 

The  development  of  a  user  oriented  interactive  computer 
program  in  solving  engineering  problems  requires  a  consider¬ 
able  amount  of  programming  work,  contributing  significantly 
to  its  complexity  and  size. 

The  computer  program,  hereafter  referred  to  as  BODPLT 
consists  of  a  main  program,  a  number  of  programmer-composed 
subroutines,  a  few  library  functions/subroutines,  and  various 
subroutines  of  the  DISSPLA  graphics  package.  The  entire 
program  is  written  in  the  FORTRAN  IV  language. 


Cascade  compensation  may  further  be  classified  into  two 
further  types;  a  high  pass  filter  usually  called  a  phase 
lead  compensator,  and  a  low  pass  or  phase  lag  compensator. 

The  selection  of  the  type  of  compensator ( s )  to  be  used 
depends  on  a  number  of  factors,  the  important  ones  being 
experience,  personal  preferences,  availability,  system  cons¬ 
traints,  etc.  Unfortunately,  there  are  no  mathematical 
techniques  to  help  in  this  selection  process.  Generally, 
one  has  to  complete  several  designs  and  then  choose  the  most 
appropriate  one. 

C.  COMPUTER  AIDED  DESIGN 

Of  all  the  classical  design  methods  available,  the  Bode 
diagram  technique  is  generally  considered  to  be  the  simplest, 
at  the  same  time  providing  the  most  insight  into  system 
performance  and  behavior.  The  Bode  design  method  may  be  used 
equally  effectively  both  with  cascade  and  feedback  compen¬ 
sation  schemes. 

The  computer  program  developed  in  this  thesis  displays 
on  the  terminal,  initially,  the  Bode  diagram  (magnitude  and 
phase)  of  the  uncompensated  open  loop  system.  The  user  may 
then  select  the  type  of  compensator  to  be  used,  and  feed  this 
information  to  the  computer.  The  display  changes,  now  showing 
the  compensated  system  Bode  Plot.  This  procedure  may  be 
repeated  iteratively,  with  the  computer  updating  the  Bode 
Plot  with  each  change  made. 


Compensators  used  in  the  design  of  feedback  control 
systems  are  generally  classified  as  shown  in  the  block 
diagrams  of  figure  1. 


COMPENSATOR 


PLANT 


o)  CASCADE  COMPENSATION 


b)  FEEOBACK  COMPENSATION 

plant  plant 


Figure  1.  Classification  of  Compensation  Structures 
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This  thesis  does  not  address  the  problem  of  converting 
from  one  form  of  representation  to  another.  It  is  assumed 
that  any  conversions  necessary  have  already  been  performed 
and  that  the  system  is  represented  by  its  open  loop  transfer 
function . 

B.  PROBLEM  FORMULATION 

Each  system  design  has  its  own  unique  characteristics, 
but  in  general  the  system  has  to  meet  some  kind  of  perfor¬ 
mance  standards.  These  performance  standards  are  generally 
provided  as  numerical  specifications.  The  first  step  in  the 
design  of  a  control  system  is  to  analyze  the  system  by  itself 
in  the  usual  feedback  loop  configuration.  This  is  usually 
referred  to  as  the  uncompensated  system. 

Analysis  of  the  uncompensated  system  almost  always  shows 
that  the  system  cannot  meet  some  or  all  of  the  given  perfor¬ 
mance  standards.  Usually,  additional  components  have  to  be 
inserted  in  the  system  for  the  purpose  of  altering  the 
performance  of  the  system.  These  components  are  called 
compensators.  Compensation  is  a  two  step  procedure,  in  which 
additional  components  (compensators)  are  inserted  to  change 
the  structure  of  the  system,  and  these  components  are  then 
adjusted  until  the  performance  characteristics  are  satisfied. 

The  theory  of  cascade  and  feedback  compensation  is 
discussed  in  detail  by  Thaler  (Ref.  2)  in  Chapters  5  and  6. 
Only  a  brief  discussion  of  the  types  of  compensators  is 
presented  here. 


EXAMPLE  4.2 


Uncompensated  Bode  Plot  of  Example 


EXRMPLE  4.2 


Table  2.  Tabular  Output  of  Example  4.2 


FRE  a 


O.iQOOOOE+OO 
0 .110  7  76  E  +00 
0.122712E  +  0J 
0  .135  935E+00 
0. 150584E+00 
0.16&810E+GU 
0.184765  E+00 
0 .2Q4697E  +00 
0  «2i.6  7  54  E  +  00 
0  .251  188E+00 
0.278256E+C0 
0 .80a  240  E  +00 
0.841454E  +  C  0 
0.378248E+00 
0 .41*  00a £+00 
0.+64158E+C0 
0  .5141 75 E+00 
0 .569  5a0.E+00 
0  •  o30  956E  +  Co 
0  .698  946E  +00 
0. 7742o3  t  +  00 
0 .857  6S5E  +  00 
Q.950117E+00 
0.1 O5250E  +  0i 
0 . 1 a6  591 E  +  0x 
0.129  X  55E+0 1 
0.143  0  72E  +  01 
0  .158+O9.E+01 
0.1 75  567E  +0i 
0.i94485E+Q1 
0  .215443  fc+Oi 
0.238658E+C1 
0  .2643  75E  +  0i 
0 . 292  664  c  +  01 
0.324421E+G1 
0  .859361  C+Oi 
0.U96106E+01 
0 .44i 0O4E+0i 
0  .4oa  526  E+01 
0 .541 166E+01 
0  .599  483  E  +0i 
0 . 664  081 t  +  Ql 
0 . 735640E  +  Cl 
0  .oi4  911  E+Oi 
0.902  723E  +  01 
0 .999  999  £  +  0i 
0  .1 1 0  7  75  £  +  02 
0.122712E  +  C2 
0  .135935E+02 
0.150583E  +  02 
0.166809E+02 
0  .1  84  7a4  E  +02 
0  .204696E  +  02 
0  .i-^6  754E  +  Gi. 
0.251 loaE+02 
0.27a255E  +  02 
0  .306239E+02 
0.3  41  454  E  + 02 
0  .878 247 E  +02 
0  .419007  £+02 
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0.  596845c  +  O2 
0. 5d*089E+G2 
0.  5  79  754E  +  02 
0.  570825E  +  02 
0. 560763E+02 
0.  551 1  06c  +  O2 
0. 541270E+02 
0. 531249E+02 
0.  52iQl2E+  02 
0.  5iG529E  +  Qi 
u.  499765E  +  02 
0.  486667E  +  02 
0.477266E+02 
0.  4e54 74E  +  02 
0.453269E+02 
0. 440698E+02 
0.  427697E  +  02 
0. 414292E+02 
0.  4Qo502c  +  02 
0. 38635oE  +  02 
0.371891E+02 
0.  357i55E  +  02 
0. 342201E  +  02 
0.  827090E  +  02 
0.  31i68lE  +  02 
0.296640E+02 
0.  28A432E  +  02 
u.2o6820E+02 
0.  251367E+02 
0.  Z3oe3iE  +  02 
0.222 164E+02 
0.206014E  +  02 
0. 194217E+02 
0. 160801E+02 
0.  167783E  +  Q*. 
U • 1551&6E+02 
0.  1+2947E  +  O2 
0.  13H07E  +  02 
0. 119624E+02 
0.  1  06‘t66E  +  02 
0. 97599 4E +01 
0. 869860E+01 
0.  765877E  +  01 
0. 66365aE  +  01 
0. 562628t+01 
0.  468006E  +  01 
0. 363853E+01 
0.  2o501 1E  +  01 
0.  166 180E  +  01 
0.  66  **65dE  +00 
-0.  88J.896E  +  00 
-0. 134252E+01 
-0.  23O847E  +  01 
-0.  84i815E  +  0i. 
-0.446056E+01 
-0.  557+SiE  +  Oi 
-0.669992E+01 
-0.  76596OE+01 
-0.  905738E  +  01 
-O. iQi96iE+02 


PrtASE 


—  0.  99  5i41 E+02 
-C.1OO5Q0E+U8 
-0.1 Oi58e  E+03 

—  0. 1027  69E+  03 
-0.104056E+U3 
-u.i03457E+08 
-0.1 069666+03 
—0.1 08594E+j3 

—  0.  1 1 08  3i  E  + 08 
-0.1 12i 74t+u3 
-0.114115E+03 
-0. 1 1613oE+u3 
—0.1 16i.26E+03 
-C.  1 2  03  5  3  E+  03 
-0.  12^.49 1E  +  03 

—  0.12  4d0oE+03 

—  0.  1 1.6663  E  +  03 
-0.126623E+Q3 

—  0.  13  04  49  E+03 

—  0.  132i  06 E  +  03 
-0. 1385 oU E+03 
-0.134785E+03 

—  0. 1 3575ot+03 

—  O.i  3o-+oi  E+03 

—  O.i3  66  84t+o8 
-0.137022E+08 
-0.  1  a66  77 £+03 

—  0.  13  645  7  E+03 
-0. i 3  57  7 3c+Oa 
-0.  U+8  53E+03 

—  0.1337 i6t  +  06 
-0.1 3239  5E+03 
-0.  130927E+U3 
-0.1 2 93 51 E+03 
-0.i2 7709 E+03 
-0.126043E+03 
-0.1 2  43  93  E+03 
-0.1c2797E  +  0-> 
-0.121289E+O3 
-O.i 1 9699E+03 
— 0.116652E+03 
-O.i 1 7569E+03 
—O.i i6e  66t+0 j 
-0.11 5955E+  03 
-0. 1i5+45E+03 
-0.113145E  +  03 

—  0.11 5o 60 E  +  03 
-0.1 15192E+03 

—  O.i  1554  7E+03 
-0. 1 16125E+03 
-0.  ii69i.7E+03 
-0.1 17953E+G3 
-0.1 1 92  ui E  +  03 
-O.it  Oo  00  E+  03 
-0 . 12 234a E+03 
-0. 1 2  4i 1 9  E+03 
-0.  126281E+03 
-O.i 26510E+03 
-0.  i3  0883E  +  03 

—  0.i333  73E  +  Uj 
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Table  2.  (Contd.) 


0  •  4o4157E*  C  2 
0  .514174E+C2 
0 . 5o9  5  79 E *02 
0 .630  955E  +  C2 
0  .696  945  E  +02 
0.774261E+02 
0 .857 692E+02 
Q.950116E-*G2 
0.1C5250E+C3 
0  .j.lo591£  +  G3 
0.129 155E+03 
0.143G72E  +  C3 
0 .158  489  E  *03 
0 . i 75  5t7  E+03 
0  .194485E+C3 
0  .215443  E*03 
0.238  6  58E  +  C3 
0  .Z64375E+C3 
0.292663E+03 
0.324421E+C3 
0  .359380E  +  C3 
Q.3961C5E+CJ 
0 .441 0C4E  +  C3 
0  .468  5^.5  E  *03 
0 .541 187E+0J 
0  .599462E  +  C3 
0.664 080 E+Qj 
0.735639E+C3 
0  .8149  j.OE  +  03 
0 . 902  721 £*C3 
0.5999S8E+03 


-0.  115780E+O2 
-0. 129043E+02 
-0.  14z75-St*u2 
-0.  156905E*  Oh 
-0. 17i484E  +  0<l 
-0. i6o46oE+02 
- 0.  <iG182cE*  Oh 
-0.217522E+02 
-0.  23352oE  +  02 
-C.  ^498C6E  +  03 
-0. 266322E+02 
-0.  28304jE  +  02 
-0. 299941E+02 
-0.31o96oE+02 
-0.334i6xE  +  Qc 
-0. 351439E+02 
-0.  ^6bo05E  +0^ 
-0.  386345E+  On 
-0.  40374 5E  +  02 
-0.  421296E  +  02 
-0. 438886E+02 
-0.456513E+O^. 
-  0. 4  74166E  +  OH 
-0.441842E*02 
-0. 509536E+02 
-0.5^7247E  +  02 
-0.  544969E  +  0^. 
-0.  56z7C3E+02 
-C.560443E+02 
-0. 5 98 190 E +02 
-0.  615944E  +  Cc 


-O.i J5950E+03 
-0. 13o56iE*03 
-0 . 141232E+03 
-0.  i.438  7l £+06 
-G.14646/E+06 
-0.1 4a*93E+03 
-0.l5l4^7E+06 
-C.153749E+03 
-0. 1 55946  E+U3 

-  C.J.58Q14E+06 
-G.l599^4£  +  0-> 
-0.1O17 Jo E+03 
-C.163392E+03 
-Q.164917E+03 
— C.16ojioE+0j 
- C. 16  759  7E  +  U3 
-0.1od765E+u3 
— C.169830E+06 
-0. 1  7  07  98E+0  ji 
-C.l  7  it>78E+Q6 
-C.1724  76E+03 
-C.173i99E+G3 
-C.  i.  7o854E+0o 
-0. 174447E+Q6 
— 0.174984E+06 

-  C.  i.  7  54o9£+06 
-0.17  54  07E+0-) 
— 0.1 766 04 £+03 
-C.176q62£+Q3 
-0.1769O6E+OJ 
-C.  l77h  76E+Gj 


3. 


and  the  tachometer  transfer  function  is  KfcS .  The  tachometer 
is  fed  back  around  all  of  the  forward  gain.  Using  Bode 
diagram  methods,  set  to  provide  =  1.3. 

b.  Solution 

Bode  Plot  for  this  system  is  shown  on  figure 
6.  The  system  has  a  phase  margin  of  about  6  degrees.  For 
Mpw  =1.3,  a  phase  margin  of  45  degrees  is  required.  A  rough 
graphical  design  on  the  uncompensated  Bode  Plot  gives 


1/H  =  12.0/S  or  H  =  0.08S 

The  Bode  Plot  for  the  compensated  system  is  on  figure  7, 
showing  a  phase  margin  of  50  degrees . 

The  close  loop  frequency  response  of  this 
example,  drawn  using  BODPLT  is  shown  on  figure  8. 
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EXAMPLE  4.3 


EXRMPLE  4.3 


FREQUENCY  [RRD.  PER  SEC  J 


^  t- 


Table  3.  Tabular  Output  of  Example  4.3 


FREw 


O.luUOCOE-Cl 
0.aa3646E-Ci 
0.129155E-C1 
0 • 146  780£- C 1 
0  .166 till)  2-Gl 
0 .189573E-C1 
0  .2i 5443E-Ci. 
0  .244  643  £-01 
0.27825©E-C1 
0  .-)16267E-Cl 
0 . 354J  81 E-Ca 
0.9O8428E-C1 
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0.  606303c  +  02 
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Table  3.  (Contd.) 
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0.3593 8GE  +  C2 
0  .'♦OS  423 E  +  0x 
0.4o4157E +  02 
0  .527 4S7E+C2 
0  .599482E+G2 
G.681250E  +  C2 
0 .7742  6.1  E+Gi 
0  .o79  919  E+02 
0 . 999  S  S5E  ♦  C2 
0  .  A  i-ic46t  *03 
0.129A54E+C3 
0.146779E+C3 
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0  .2A544JE+C3 
0  .244 643 E *03 
0.278255E+C3 
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0 .408 422E *C3 
0.464157  E *03 
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0  .681<c90E+03 
0 . 77426G£*C3 
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0.999997E+C3 
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-0.  644466E  +  0a 
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-0 • 1666 66 E +03 
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-C. 1 7 1179E+03 
-0.1 7 2224 E +03 
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—  C • 1 7  39  64t+0 J 
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-0.1 77189E+03 
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—  Q.  x  78o9 

-0. 1 7oo50fc+03 
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-C.179309E+OJ 
-0.  17  939AE+U6 
— 0.a79465E+03 


Response  of  Example 


V. 


CONCLUSION  AND  RECOMMENDATIONS 


A.  CONCLUSION 

The  objective  of  this  thesis  was  to  develop  an  inter¬ 
active  user  oriented  computer  program  which  would  aid  in 
solving  control  engineering  problems  using  the  Bode  method 
of  design.  The  presented  program  proves  that  frequency 
domain  design  of  control  systems  using  the  digital  computer 
as  an  aid  is  not  only  feasible  but  highly  desirable.  The 
results  obtained  are  readily  interpretable  and  provide  good 
and  meaningful  insight  into  the  problem. 

The  results  obtained  during  the  investigation  of  the 
program  performance  show  that  a  complicated  but  well  formu¬ 
lated  problem  can  be  solved  with  ease,  and  the  solution  is 
obtained  with  speed,  accuracy  and  precision. 

The  entire  program  is  less  than  1000  lines,  with  a  total 
of  9  subroutines.  Every  effort  has  been  made  to  keep  the 
program  simple  yet  unambigous,  so  that  the  user  has  to 
invest  very  little  time  learning  how  to  use  it.  Effort  has 
also  been  made  to  minimize  the  use  of  the  computer  CPU  time. 
However,  expenditure  of  CPU  time  is  to  a  large  measure 
dependent  upon: 

(1)  the  order  of  the  system. 

(2)  the  number  of  iterations  used  in  reaching  a 


(3)  the  type  and  order  of  the  compensator (s )  used. 


B.  RECOMMENDATIONS 

The  program  as  presented  in  this  thesis  seems  to  be  able 
to  adequately  satisfy  the  usual  needs  in  a  control  system 
design  problem.  A  number  of  useful  extensions  to  the  work 
developed  in  this  thesis  can  be  carried  out.  These  are 
briefly  discussed. 

1 .  Curve  Fitting 

Although  not  specifically  worked  on  in  this  thesis, 
the  program  can  be  used  quite  effectively  for  curve  fitting 
purposes.  This  was  demonstrated  in  the  initial  stages  of  the 
development  of  this  program.  The  procedure  is,  by  its  very 
nature,  iterative  and  therefore  time  consuming  and  cumbersome. 
However,  the  entire  algorithm  can  be  automated  using  a 
minimization  subroutine,  with  the  program  outputting  a 
polynomial  to  fit  a  given  curve  over  a  specified  range  of 
the  independent  variable. 

2 .  Computer  Selection  of  Compensators 

A  suitable  minimization  routine  such  as  Box  PLX 
can  be  incorporated  into  the  program  which  could  select  the 
best  possible  location  of  poles  and  zeroes  to  meet  given 
performance  specifications.  This  would  automate  the  entire 
Bode  design  procedure,  the  user  then  having  the  option  of 
only  specifying  the  type  of  compensation,  i.e.,  cascade  or 
feedback.  It  may  however  be  pointed  out  that  minimization 
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routines  by  the  very  nature  of  their  operation  are  very  time 
consuming  and  wasteful  of  CPU  resources. 

3 .  Root  Locus 

Most  of  the  subroutines  developed  in  this  thesis 
are  very  general  in  nature  and  can  be  adapted  quite  easily  to 
develop  a  similar  interactive  program  for  Root  Locus  plots. 

4 .  Integrated  Control  System  Design 

No  meaningful  design  of  control  systems  is  complete 
without  finally  analyzing  its  time  domain  performance.  It  is 
therefore  considered  highly  desirable  to  incorporate  into 
this  program,  an  interactive  Root  Locus  design  procedure  and 
then  a  time  domain  analysis  of  the  compensated  system.  The 
entire  package  would  then  be  an  excellent  teaching  aid  for 
control  system  design. 
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